Elasticity of cross-linked semiflexible biopolymers under tension 
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Aiming at the mechanical properties of cross-linked biopolymers, we set up and analyze a model 
of two weakly bending wormlike chains subjected to a tensile force, with regularly spaced inter- 
chain bonds (cross-links) represented by harmonic springs. Within this model, we compute the 
force-extension curve and the differential stiffness exactly and discuss several limiting cases. Cross- 
links effectively stiffen the chain pair by reducing thermal fluctuations transverse to the force and 
alignment direction. The extra alignment due to cross-links increases both with growing number 
and with growing strength of the cross-links, and is most prominent for small force /. For large /, 
the additional, cross-link-induced extension is subdominant except for the case of linking the chains 
rigidly and continuously along their contour. In this combined limit, we recover asymptotically the 
elasticity of a weakly bending wormlike chain without constraints, stiffened by a factor four. The 
increase in differential stiffness can be as large as 100% for small / or large numbers of cross-links. 

PACS numbers: 87.10.Pq, 87.15.La, 82.37.Rs, 36.20.Ey 



I. INTRODUCTION 

Many important biopolymers, such as DNA, the cy- 
toskcletal filaments (filamentous (F-)actin, microtubules, 
intermediate filaments), as well as collagen in the ex- 
tracellular matrix are fluctuating macromolecules with 
a bending stiffness intermediate between that of a ran- 
dom coil (Gaussian chain) and a rigid rod. Polymers 
whose elastic behavior is dominated by their bending 
rigidity are known as semiflexible. Numerous experi- 
ments probing their elasticity have become available [TJ [2] 
with the advances in single-molecule manipulation, par- 
ticularly for DNA. Intriguing and qualitatively novel me- 
chanical behavior arises if semiflexible polymers are pair- 
wise permanently cross-linked. The elasticity of cross- 
linked biopolymers is widely studied experimentally via 
force-extension measurements. In this article, we study 
analytically the force-extension relation of an irreversibly 
cross-linked pair of semiflexible polymers within a meso- 
scopic theoretical model. 

Ubiquitous as extracellular mechanical support is the 
connective-tissue protein collagen, whose fibrils achieve 
their strength via covalent intermolecular cross-links be- 
tween triple-helical molecules [31 H]. Atomic force mi- 
croscopy [S]-(7] has been used to analyze single collagen 
fibrils, which themselves consist of many microfibrils and 
hence can be modeled as anisotropic networks of irre- 
versibly cross-linked semiflexible polymers 8 . Cell shape 
and stability is provided by the actin cytoskeleton, a net- 
work of cross-linked F-actin ranging in morphology from 
a dilute mesh to bundles of parallel filaments [9 . The 
special elastic properties due to cross-linking, closely re- 
lated to biological function, thus have become a subject 
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of increasing interest and vigorous research activity [TO] ■ 
Yet, theoretical understanding is incomplete and expla- 
nations based on semi-microscopic descriptions are rare. 

Crucial experimental results such as the strong stretch- 
ing of double-stranded DNA [TT] have been successfully 
explained by the theoretical force-extension relation for 
a weakly bending wormlike chain |12j . The wormlikc 
chain (WLC) [LMT5] maps the conformations of an in- 
extensible semiflexible polymer to one-dimensional paths 
whose statistical weight penalizes curvature and is de- 
termined by two length scales only: the total contour 
length L and the directional correlation or persistence 
length L p , proportional to the bending rigidity k. The 
weakly bending approximation of a WLC |12) simplifies 
analytical treatment by assuming that the tangent vec- 
tor at any arc-length position and the end-to-end vector 
make a small angle. This approximation applies to poly- 
mers with a large persistence length L p (compared to 
L) or subjected to strong stretching. Inhomogeneities or 
inter-molecular interactions have been the subject of sev- 
eral modifications and extensions of the weakly bending 
WLC: In [IB], the force-extension relation of two par- 
allel aligned, weakly bending WLCs with a single irre- 
versible cross-link and the elasticity of an anisotropic 
network of aligned chains have been analyzed. In the 
wormlike bundle model, an arbitrary number of regularly 
arranged parallel filaments is effectively cross-linked by 
a coarse-grained, continuous interaction [TJ]. The effect 
of spontaneous polymer curvature has been studied in 
[151 IT9] . Weak extensibility of semiflexible polymers at 
strong stretching has been addressed with a combination 
of a WLC and a Gaussian chain, the semiflexible har- 
monic chain (SHC), in [20j . 

In this work, we consider the elasticity of two iden- 
tical weakly bending WLCs connected by an arbitrary 
number of cross-links regularly spaced along the poly- 
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mer contour. The cross-links are represented by entropic 
harmonic springs, which allow for a finite extent of the 
inter-polymer distance at the cross-link sites. In the case 
of infinite spring strength, we obtain the limit of hard 
cross-links (strong topological constraints) . By introduc- 
ing infinitely many cross-links at fixed contour length, we 
can also model a continuous cross-linking or an attrac- 
tive inter-molecular interaction. The ladder structure of 
our system is reminiscent of the base-pair sequence of 
double-stranded DNA, but for the reversible hydrogen 
bonding. This obvious modification of our model to re- 
versible and/or sectional cross- linking may prove versa- 
tile for future studies of, e.g., the denaturation of DNA. 

The paper is organized as follows: In Section [TTJ we in- 
troduce the model and the observable. In Section \LU\ we 
present the main steps in calculating the canonical parti- 
tion function, from which all equilibrium quantities can 
be derived. Details of this calculation are given in Ap- 
pendices [A] and [B] In Section [TV] we present the central 
result which is the force-extension relation. After pre- 
senting the general result, we particularly focus on the 
limit of hard cross-links, of continuous cross-linking, the 
linear elasticity for small forces, and the strong stretching 
limit. For the limit of continuous cross-linking, the gen- 
eral result and a short comparative discussion are given 
in Appendix [Cj We conclude and discuss further exten- 
sions of this work in Section fVl 



II. MODEL 

We consider two weakly bending, semiflexible chains 
of equal bending rigidity and contour length, aligned 
parallel along a preferential direction x and cross-linked 
at equidistant arc-length positions specified below. The 
chain configurations are described by paths J"i(s), »"2(s), 
with s € [0, L] the arc-length parameter, and tangential 
vectors tj(s) = d s rj := drj/ds. Our setup, taking into 
account space dimension d — 2 only, is sketched in Fig. [T] 
for 3 cross-links. As indicated, we assume hinged-hinged 



» cross-links 
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FIG. 1. (Color online) Stretched, weakly bending (see main 
text for explanation) chain pair connected by 3 cross-links. 

boundary conditions, implying confinement of the verti- 
cal positions (here, to y = 0) and vanishing curvature at 
the ends. These boundary conditions are motivated by 
the following situation: Experimentally, a tensile force 
can be applied via optical or magnetic tweezers that con- 
trol the position of beads attached to the polymers' ends, 
cf., e.g., [21]. Optical tweezers usually restrict the bead's 



transverse motion, but not the rotation, so that no mo- 
ments are exerted at the ends. Additionally, we assume 
xi(0) = a; 2 (0), in order to exclude an overall x shift be- 
tween the chains. 

The effective bending potential of semiflexible chains 
(without taking into account torsions |21j ) is 

n hcnd = ^ / d s ^|d s t,( s )| 2 , (i) 

with the bending rigidity k, related to the persistence 
length L p in d = 2 via 

k = ^k B TL p , (2) 

and with the local inextensibility constraint 

\t 3 (s)\ = l, se[0,L], j = 1,2. (3) 

In order to account for inextensibility in a mathemat- 
ically tractable way, we consider the weakly bending ap- 
proximation: The chains' tangents preferentially align 
with a given direction, here x. A stretching force of 
strength /, acting on both ends of the chains (cf. Fig.JT]), 
is described by the potential 

2 

Stretch = -fe x ■ (rj (L) - Tj (0)) (4) 

J'=l 

2 ,L 

= -/Y" / dsd s Xj(s). 

For sufficiently large stretching forces or bending rigidi- 
ties, the tilt of the tangent vector away from the x axis is 
small, so that the condition Eq. ^ reads approximately 

d sXj (s) = \-\ (d s y 3 (s)) 2 + O ([dsVMf) ■ (5) 

Inserting this expansion into Eqs. (JIJ and Q and dis- 
carding all but quadratic terms in derivatives of y, we 
arrive at the weakly bending approximations of the bend- 
ing and stretching potentials. 

Cross-links between the two chains are introduced at 
N — 1 sites regularly spaced along the contours, 

s b = b ^, b=l,2,...,N-l, (6) 

dividing the contour length L into N sections, cf. Fig.[T] 
Explicitly, we model cross-links as entropic, harmonic 
springs of strength g = 2k B Tja} c , where a 2 c is the 
temperature-independent squared equilibrium length of 
one cross-link. 

Finally, the total effective Hamiltonian is 




Wo 

N-l 



+ I^2(yi(s»)-y2(s b )) 2 , (7) 

6=1 
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where "Ho is the Hamiltonian of the system without cross- 
linking, and the last term oc k^T describes the entropic 
cross-links in weakly bending approximation. 

Starting from the concept of harmonic cross-links at 
discrete sites Sb = bL/N, we will also consider the limit of 
continuous cross-linking, achieved by taking N — > oo and 
As := L/N — > 0. In this case, the strength g of a single 
cross-link has to go to zero, such that the total strength 
g '■= Ng remains finite. Replacing Yl^S^ — > j; m 
the cross-link part of Eq. ^ gives a continuous, harmonic 
inter-chain attraction of strength g/L, 
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ds (yi(s) - y 2 (s)Y 



(8) 



It is our aim to study the effect of cross-links on the 
chain elasticity and hence to compute the force-extension 
relation exactly for an arbitrary number of irreversible 
cross-links. Thus, the relevant quantity is the average 
end-to-end extension of one chain in force direction x, 



(x) := (x(L)-x(0)) n = L \ ^J\s(d s y) 2 ^ , 



(9) 



where y/n, denotes the canonical average with the Hamil- 
tonian of Eq. ([7]) . 

The force-extension relation of one weakly bending 
WLC without cross-links (Hamiltonian Ho), first ad- 
dressed by Marko and Siggia [T5], is 



( x )n n =1 L f coth^ 1 \ 

L 2L p [ y/Jr fr J 

in terms of the dimensionless variable 



f r := fL 2 /K, 



(10) 



(11) 



which is the ratio of stretching energy, /L, and bending 
energy, k/L. The force-extension curve displays a lin- 
ear regime for small forces / and in the limit of strong 
stretching approaches the maximal end-to-end extension 
L with a characteristic saturation oc J -1 / 2 . 



III. PARTITION FUNCTION 

In this section, we detail the calculation of the canoni- 
cal partition function, Z — J V[y(s)} e — ^foWJ (the con- 
figurational integral for both chains denoted by T>[y(s)]), 
which provides access to all equilibrium observables. For 
the purpose of this work, the end-to-end extension de- 
fined in Eq. ^ is obtained from InZ or the free energy 
[2"2"] by differentiation with respect to the force /: 



(x) 



k B Td\nZ 

~2 dT' 



The first step is to expand the chain configurations 
yj(s) in appropriate eigenfunctions. As mentioned above, 



we impose hinged-hinged boundary conditions, which for 
our system translate into [y'J (s) := d^y(s)] 

yj (0)= yj (L)=0 y'JiO) = y'J(L) = 0, 3 =1,2. (13) 

According to these boundary conditions, our Fourier- 
series ansatz is 



M 



2/i («) 



^ A m sm(q m s), 

m— 1 
M 

2/2 (s) = X! B m sm(q m s), 



with wave numbers 

TO7T 



L ' 



771=1 



m the mode number, 



(14) 



(15) 



and M the largest undulation mode considered within 
our continuum model (roughly, the wave-length reso- 
lution is bounded by molecular distances). With this 
ansatz, the Hamiltonian H can be written as a quadratic 
form in the coefficient vector 



T := (A U B U A 2 ,B 2 ,...) T . 
Omitting the constant —2fL, 



(16) 



A I 



'=1 



where C, due to T-Lq of the uncross-linked system, is a 
diagonal matrix ((g) denotes the Kronecker product), 



C = diag(ci,c 2 , 

Cm := j(f + Hq 2 mYl 



1 
1 



2 

in " 



(18) 



and UU T the matrix due to the cross-link Hamiltonian, 
cf. the second line of Eq. Q, or Eq. 

The partition function follows as a generalized Gaus- 
sian integral over the mode coefficients normalized by L, 



Z = Jv[T] e-wm, 



(19) 



M 



V[T] := J] d 



A, 



B„ 



Returning to the end-to-end x-extension introduced in 
Eqs. Q and ( 12 ), we wish to focus primarily on the cross- 
link contribution, i.e., 



(12) (Ax)i=(x(L)-x(0)) n -(x(L)-x(0)) no , (20) 



since the extension (x(L) — x(0)^^ due to thermal fluc- 
tuations of uncross-linked weakly bending WLCs only is 
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known . To that end, we write the partition function 
as Z = Z Te \Za, where Zq is the partition function of the 
uncross-linked system and address the relative partition 
function [cf. Eq. p7|], 



-Tel 



JV[T]e-^ 
~~ JV[T] e-^«o[r] 

(dettl + C- 1 ^))- 172 . 



(21) 



This yields the cross-link-induced extra displacement 

k B T 91nZ rel 



(Ax) 



Of 



(22) 



A. Finite number of cross-links 

First, we address a finite number N — 1 of equidistant 
harmonic cross-links, for which the cross-link Hamilto- 
nian is quadratic, but not diagonal in the modes. The 
matrix UU T is a sum of N — 1 projectors, 



N-l 

^ T = § E K ® O 

6=1 



1 -1 
-1 1 



(23) 



(sin(gis b ),sin(g 2 s b ), . . .). 



Using the identity det exp A = exp tr A to expand the 
determinant in Eq. (21), Z le \ is accessible via traces of 



powers of the matrix C UU 



Z rcl = exp 



I k=l 



(_l)fe+l 



T\k 



(c-'uu 1 ) 



In the trace of a power k of C 1 UU T , the 2x2 matrices in 



the Kronecker products, Eqs. (18) and (23), merely pro- 
duce a factor 2 fc , which can be computed separately be- 
fore performing the trace operation over the mode indices 
rrij. Hence, we are left with handling the non-diagonal 
mode-index structure of the projector sum Eq. (23 1. The 
corresponding matrix of rank (N — 1), 



JV-l 



T 
b ■ 



(25a) 



6=1 



with mode indices mi , mi has entries 



jV-l 



P„ 



N 

5., 



sin I sm 



6=1 

mi—m 2 ,2ZN 



N 



N 



'mi+m 2 ,2ZA r ! 



(25b) 



where Z denotes the set of integers, such that P is a 
sparse matrix of block form: In each quadratic block of 
dimension 2N, nonzero entries appear on the diagonal 
(+1) and on one anti-diagonal (—1) only. For the explicit 
form of P, see Appendix |A} The rows/columns display 



the occupation structure of the N — 1 eigenvectors (la- 
beled by / in the following) in the Fourier basis, each with 
nonzero amplitudes only for a subset of modes, which at 
all cross-link sites are pairwise in-phase or phase-shifted 
by 7r. Modes indexed by multiples of N have nodes at all 
cross-link sites, thus do not contribute to the cross-link 
energy, and constitute the kernel of P. 

Due to this special form of the matrix C~ 1 UU T , we 
are able to derive a closed expression for the trace of any 
power of C~ X UU T (cf. Appendix |b|) , viz. 



f) "tr(<rW)* 




(26) 



C 2{ti-l)N+l + C 2nN-l 



Here, due to the fast decay wit h m ode number of the 
inverse coefficients c" 1 from Eq. ( 18 ) - basically inverse 
elastic constants for the undulation modes - we have ex- 
tended the summation over modes to a series. Combin- 



ing Eqs. (24) and (26), we hnd that the relative partition 
function Z Ie \ factorizes into N — 1 different "eigenvector" 
factors, or equivalently, 



In Z 



rcl 



(27) 



w-i 



2 ^ 

l=i 



\n< 1 



gN 
~2~ 



OO 

H ( C 2(M-1 



-DAM 



-2fiN-l 



JV-l 

(=1 



(24) By inserting the into the series, we obtain for the 



factors of the partition function, Eq. (27) 



Zr = {l 



sinh 5 f 



'1/2 



(28) 



cosh Sf — cos <f>i ' 
Herein, we employ the dimensionless variable 



(29) 



which is the ratio of two lengths: The arc-length spacing 
L/N between cross- links and the directional 'memory' 
length WkJJ of the stretched WLC, or the penetration 
depth of boundary conditions [TJ]. The dependence of 
the partition function on the bending rigidity k is via this 
ratio only. Additionally we define the phases specific to 
the N — 1 eigenvectors, 

<h ■■= ^- (30) 

We note that the Gaussian statistical weights and the 
regular cross-link spacing simplify enormously, if not en- 
able at all, analytical calculations. 
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In fact, Eqs. (27) and (28) are a central result of our 



paper, yielding the exact free energy of two cross-linked, 
weakly bending WLCs as F = F + AF, where F is the 
free energy of the chain pair without cross- links, 



M 



T 2 r 



F = -k B T In Z = k B T > In — , (31) 

m— 1 



(leading to the force-extension relation Eq. (10)), and 
the free-energy increment AF due to cross-links is the 
sum of the Fi = -fc B TlnZ ; from Eq. Via Sf, this 

free energy depends on the dimensionless energy ratio 
f r = N 2 tfj introduced in Eq. (jTTJ) . As already mentioned, 
use of the weakly bending approximation requires that 
either the work done by the external force / or the bend- 
ing energy is large compared to the thermal energy, in 
order to restrict the transverse fluctuations to be small. 
The dependence on the free parameter f r will be fur- 
ther discussed for the cross-link contribution to the force- 
extension relation. In addition, the free energy depends 
on the ratio of cross-link energy to stretching energy, 



_ gL 

69 ' Nf 

and, of course, on the number of cross-links, N — 1. 



(32) 



B. Continuous cross-linking via harmonic 
inter-chain attraction 



Here, we sketch the derivation for an infinite number 
of regularly spaced cross-links, N —> oo, at finite total 
strength g := Ng, and for finite contour length L. With 
the continuous cross-link Hamiltonian Eq. (|8| in normal- 
mode representation, 

M 

U{C) = il E (A mi -B mi )(A m2 -B m2 ) (33) 
x / ds sm(q mi s)sin(q m2 s) 



JO 
. M 



m—l 



the total Hamiltonian is diagonal with respect to the 
mode indices. Thus, we find for the excess free energy due 
to the inter-chain attraction a closed expression, again 
extending the sum over the modes to a series, 



AF (c) = -k B T In Z^] 



^ V in r I " 



(34) 



in agreement with performing the limit TV — ¥ oo at finite 



IV. FORCE-EXTENSION RELATION 



Using Eqs. (27) and (28), straightforward yet tedious 
differentiation with respect to / yields the force-extension 
relation 



(Ax) 



k B Tg sri ni{5 f ) 
8V/2 ^ di(8 f) e g y 



with numerator 



ni(Sf) =2 (coshSf — cos (j>i) 



(35) 



(36a) 



- (1 - cos 
and denominator 



j sinh Sf 1 — cos 0; cosh S j 
Sf cosh Sf — cos (f>i 



di (5f,t 



(1 — cos0;) (cosh Sf — COS tfil) 



cosh Sf — cos (f>i — (1 — cos 4>i) 



(36b) 



sinh S 



V 



in terms of the length ratio Sf from Eq. (29), the ratio 
e g from Eq. ( 32 ) , and the phases <fii from Eq. ( 30 1 [23 . 



Since a direct interpretation of the expressions in 
Eqs. (|35) and @ is difficult, in Fig. [2] we show the ca - 
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FIG. 2. (Color online) Force-extension relation at a relative 
persistence length L p /L — 10: for different numbers of cross- 
links at finite cross-link strength g a — 50 (top) and for differ- 
ent cross-link strengths at 10 cross- links (bottom). 
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culated force-extension relation for several numbers and 
strengths of cross-links, as a function of the dimensionless 

The dimensionless parameter 
,h is 



force variable f r , Eq. (11 



for the cross-link strengt 



9a 



knT 



L 2 



(37) 



which, by virtue of our entropic-spring model for the 
cross-links, can be expressed as the squared ratio of the 
WLC contour length and one cross-link's length at rest. 

First, in the upper part of Fig. [2j the force-extension 
curve is plotted for different numbers of cross-links, at 
constant strength of a single cross-link. The overall form 
of the saturation curve is reminiscent of an unconstrained 
weakly bending WLC [12]. Evidently, the general effect 
of cross-linking is to increase the extension in force direc- 
tion relative to an uncross-linked weakly bending chain, 
because cross-links effectively suppress thermal fluctua- 
tions perpendicular to the aligning force. The growth 
of the extra alignment with the number of cross-links 
is nonlinear, the increase relative to a chain pair with 
less cross-links being largest for a few cross-links, and for 
weak stretching. The limit of continuous cross-linking, 



cf. Sec. (IIIB), is discussed in Sec. IV B 



Enforcing a smaller and smaller cross-link length (in- 
creasing the cross-link strength) enhances the alignment 
or effective stiffness, too, as visible in the lower part of 
Fig. [2j in which the cross-link strength is varied at a 
fixed number of cross-links. The limit of strong topologi- 
cal constraints at the cross-link sites (hard or inextensible 
cross-links) is presented in Sec. IV A| 

Cross-links are most effective in suppressing transverse 
fluctuations and aligning the chain pair at small reduced 
stretching forces, at which the directional memory length 
yj 'k/ f is still large compared to the cross-link spacing 
L/N. For these relatively weak pulling forces, there is a 
regime of linear elasticity for all numbers and strengths 
of cross-links, taken a closer look upon in Sec. |IVC| 
For increasing force, the incremental extension due to 
cross-links decreases, since at strong stretching, the dom- 
inant contribution to the saturating extension arises from 
"pulling out" the remaining length reserves stored in 
thermal undulations. The asymptotic decay of the cross- 
link contribution with force is computed in Sec. |IVD| 

In Fig. [2j we have chosen a ratio of persistence to con- 
tour length Lp/L — 10 sufficiently large as to give for all 
f r relative extensions close to 1, in order to explore the 
entire range of stretching forces and yet keep the weakly 
bending approximation. A ratio L p /L of this order would 
apply to long microtubules [21] ■ The persistence length 
of actin is about 15 [im [22[[25], thus for typical lengths 
of actin filaments in solution, the ratio L p /L is of or- 
der 1. For smaller ratios, e.g., L p /L ~ 0.1 for type I 
collagen fibrils [26l [27], or L p /L ~ 0.01 for 10 fim of 
double-stranded DNA [TT], our predictions are reason- 
able at strong stretching only. 



A. Limit of hard cross-links 

In the limit of infinite cross-link strength or vanishing 
ratio of cross-link to contour length, a c /L — > 0, we have 



with ni(5f) from Eq. (36a 



(AsW) = faT^! ni (S f ) 
and 



dj (6f) — cosh Sf — cos 4>i — (1 — cos <pi) 



sinh 5 f 



(38) 



(39) 



The corresponding force-extension curves are shown in 
Fig. § for more flexible weakly bending chains (with 



hard X-links, persist, length L p = L 




FIG. 3. (Color online) Force-extension relation for different 



numbers of hard cross-links and L p = L. 



te rels 



and 



L p = L) than the rather rod-like chains in F 
[I] Since hard cross-links completely eliminate relative 
motion of the filaments transverse to stretching at the 
cross-linking sites, the relative alignment effect due to 
cross-links is seen to be stronger. 



B. Continuous cross-linking 

In this Section, we discuss the limit of continuous cross- 
N 



linking from Sec. IIIB N — >• oo at total strength g := 
Ng. The incremental extension of one chain due to con 



tinuous cross- linking is computed from AF^ C \ Eq. (34 1, 



by differentiation, cf. the general expression for all values 
of g and further remarks in Appendix [C] In the case of 
continuous and rigid cross-linking [a c — > 0), the force- 
extension relation is 



L 



= 1 



L ( coth 



1 

Tr 



(40) 



Comparing this result to Eq. (10), we observe that the 
squared thermal y-fluctuations, cf. Eq. Q, are reduced 
by 1/2 relative to the uncross-linked case. This does, 
however, not imply that the two chains attached to each 
other rigidly can be treated as unconstrained weakly 
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bending WLCs with just one effective persistence length 
or bending stiffness n c g, since f r itself depends on the 
bending stiffness k. The different apparent persistence 
lengths in the force regimes of linear elasticity and of 
strong stretching are discussed in the next two Sections. 
The extension for rigid, continuous cross-linking for ar- 
bitrary force, Eq. ( 40 ) , is shown as the topmost curve in 
Fig. |3j corresponding in Fig. [2] to the asymptotic case of 
both infinite number and strength of cross-links. 



C. Force- free extension and linear response regime 

Knowing the exact extension curve for all values of the 
force / allows us to address the equilibrium extension 
in the limit / — > and the linear elasticity for small / 
- assuming a large persistence length L p , so that weakly 
bending holds. This linear response or weak-perturbation 
regime may be the best accessible for stretching exper- 
iments on sensitive biopolymers. Moreover, the force- 
extension curves computed within our model suggest that 
the chain pair's extension at moderate or zero force, 
cf. Fig. [2j is most indicative of the degree of cross- 
linking. Without cross-links, the equilibrium extension 
of a weakly bending WLC parallel to alignment for / —> 
is, following Eq. (10 1, 



L 



lim 

/->o 



L 



L 
6L 



(41) 



v 



cf. [21] . The deviation from the maximal extension is 
inversely proportional to the persistence length. Cross- 
links increase the equilibrium extension according to 



L 



(42) 



N-l 



x 2 + 13xi + W] /(1-xi) 
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tt 6(l-*0 + 



with xi := cos (j>i. For large L p and finite g a , this expres- 
sion is O (L/Lp) 2 and hence small. For hard cross-links, 
Eq. (42) is linear in L/L p , so that an effective persis- 
tence length L PjG ff > L p of the cross-linked chains can be 
defined, viz., 



= 1 



N-l o 



13xi + 16 



J p,cS 



20N 2 ^ (l-xi)(2 + x l Y 



(43) 



In the limit of continuous, rigid cross-linking discussed in 
Sec. |IVB[ the increase in equilibrium extension is 



<A4 c ' h) ) 



L 



12L V 



(44) 



Upon comparison with Eq. (41), we thus find the zero- 



force extension of one weakly bending WLC with twice 
the original persistence length or n ff = 2n. 
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FIG. 4. (Color online) Linear elastic constant as a function 
of the number of cross-links, N — 1, for L p /L — 10. 



In Fig. |4j we show the dependence of the linear elastic 
constant, computed in dimensionless form as 



E Q := dj 



(45) 



/r-=0 



on the number of cross-links. The elastic constant of an 
uncross-linked weakly bending WLC is given by Eq — 
90L p /L. For extensible cross-links, a large increase with 
cross-link number up to N « 10 is followed by a satu- 
ration to twice the elastic constant of an uncross-linked 
weakly bending WLC. For hard cross-links, the increase 
in the elastic constant caused by introducing only a few 
cross-links is even more drastic, and the curve approaches 
a step function. 



D. Strong stretching limit 

Here, we consider the limit of strong stretching, i.e., 
f r ^> 1, or for finite N, \ffr/N ^ 1, which means that 



the directional memory length introduced after Eq. ( 29 ) 
is much smaller than the cross-link spacing, 



k L 
1 <<: N- 



(46) 



At finite cross-link strength and for a finite number of 
cross-links, the dependence on the individual eigenvector 
phases remains in the limit of strong stretching, yet the 
asymptotic scaling with f r is the same for all summands, 



(Ax) 
L 

= g a L 2 
NLj 



(47) 
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The same asymptotic decay ensues for continuous cross- 
linking at finite total strength g a := 2NL 2 /a 2 cl cf. 
Eq. (|C1|), but with a simpler coefficient, 



3L2 



fr + O 



(/; 



5/2 



(48) 



In the presence of a finite number of hard cross-links, 
strong stretching asymptotically results in an extension 
increment, which is independent of the individual eigen- 
vector phases and decays proportional to /~ , 



hard X-links, persist, length L D = 10L 



(AsPO) _ (N-1)L j 
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2L r . 



Of, 



r-3/2 



(49) 



Hence the asymptotic decay of the extra alignment with 
force is markedly slower than for extensible cross-links. 
In all cases mentioned so far, the impact of cross-linking 
diminishes fast for strong stretching, and the extension 
curve displays the saturation cx f^ 1 ^ 2 of the uncross- 
linked chain's extension to the contour length. 

Obviously, the asymptotic behavior computed for hard 



cross- links in Eq. ( 49 1 cannot apply in the limit N — > oo 
of continuous, rigid cross-linking, due to the diverging 
prefactor oc (N — 1). Indeed, in this case, the asymptotic 
decay of the inter-chain contribution to the extension is 
even slower, viz. 



L 



L 

4Lr. 



f, 



-1/2 



O (fr 1 ) 



(50) 



This is the same scaling as the saturation of the uncross- 
linked chain's extension, hence the stabilizing effect of 
continuous, rigid cross-linking is manifest even for large 
stretching forces. Moreover, in the strong stretching 
limit, a continuously and rigidly linked chain behaves ef- 
fectively like a weakly bending WLC with fourfold origi- 
nal persistence length or K c g = 4k. 

Of course, the WLC picture is oversimplified at very 
strong stretching, at which inextensibility is clearly vio- 
lated for many semiflexible biopolymers [H [25] . 



E. Differential stiffness 

A quantity of interest related to the x extension is the 
differential stiffness: In the corresponding experiment for 
our system, both fibers are pre-stressed by longitudinal 
stretching, subsequently, the strain response to a small 
change in the applied stress is measured. From the force- 
extension relation, we can readily compute the differen- 
tial stiffness as the quotient of force and extension in- 
crement at a given pre-stretching force. More precisely, 
we consider the dimensionless differential stiffness as a 
function of the dimensionless force f r , 



E(fr) ■= d fr 



(51) 



generalizing the elastic constant discussed in Sec. |IV C| 

In an effort to highlight the effective stiffening due to 
cross-links, we show in Fig. [5] the increase in differential 
stiffness relative to uncross-linked weakly bending WLCs 
for the case of hard cross-links. Again, for a few cross- 
links, the differential stiffness is particularly enhanced 



c 

to 

TD 

l/i 
CO 



100 




80 




60 


" X4inks: 
1 

- 2 


40 


20 


4 

- 10 





100 



10 



20 



30 



40 



50 



force t = fL /k 



FIG. 5. (Color online) Increase in differential stiffness relative 
to a weakly bending WLC without cross-links. 



at weak stretching. Already for a single, hard cross-link, 
the linear elastic constant is increased by more than 80 % 
compared to the uncross-linked case. Upon approaching 
the limit of rigid, continuous cross-linking, cf. Eq. (40 1, 



the differential stiffness is increased by a factor two for 
all values of f r . 



V. DISCUSSION AND OUTLOOK 

Within a transparent mesoscopic model of cross-linked 
polymers, the elasticity of two irreversibly cross-linked 
WLCs subjected to a tensile force has been studied in 
the weakly bending approximation. The validity of the 
latter is granted by assuming either a large tensile force 
or a large bending rigidity. For an arbitrary number of 
cross-links with given strength, we have calculated the 
free energy and, derived thereof, the force-extension re- 
lation exactly. Both with increasing number N and with 
increasing strength g of the cross-links, the effective stiff- 
ness of the chain pair increases, since cross-links stabilize 
the chains against thermal undulations. Particularly for 
weak stretching, the enhancement in alignment is con- 
siderable, such that in corresponding weak-perturbation 
experiments on biopolymers, the increase in the linear 
elastic constant may be a useful indicator of (partial) 
cross-linking. As expected, the effect is most pronounced 
for hard cross-links. In the limit of strong pulling forces, 
the additional extension (Ax) due to cross-linking de- 
creases, and the elasticity of an uncross-linked WLC [12] 
dominates. However, the asymptotic behavior for large 
stretching forces is different for hard and extensible cross- 
links, as well as for discrete and continuous cross-linking, 
and is summarized in Table [I] For extensible cross-links, 
the cross-link contribution decays as / _2 , for a finite 
number of hard cross-links, as / . A slower decay is 
found in the limit of both cross-link number N — > oo 
and cross- link strength g — > oo, in which the two chains 
are linked continuously and rigidly along their contour: 
For strong stretching, the asymptotic form of the force- 



9 



extension relation reflects the behavior of one uncross- 
linked weakly bending chain with effective persistence 
length 4L p and with the known J -1 / 2 scaling. 

TABLE I. Exponents of the asymptotic scaling of (Ax) with 
force f r = fL 2 /n in the strong stretching limit, f r 3> 1. 

number of cross-links 
finite infinite 
extensible —2 —2 

hard —1 -1 /2 



method two parallel aligned WLCs whose arc-length is 
sectioned into cross-linked and disconnected parts. In 
the set-up we have considered here, all ingredients, viz., 
cross-linking, bending stiffness, and longitudinal forcing, 
act to decrease transverse fluctuations of the chains. An 
unzipping transition could presumably be studied in an 
altered situation, e.g., one, in which the cross-linked fil- 
aments are teared apart at one end. More complicated 
refinement of our model might account for twist and over- 
stretching, effects shown to be essential for the elasticity 
of DNA [37]. 



From the exact extension for all stretching forces, we 
have computed another experimentally relevant observ- 
able, viz., the differential stiffness of the (pre-stretched) 
cross-linked chains. Even a small number of cross-links 
enhances the differential stiffness dramatically. Again, 
the impact is largest for small stretching forces, which 
can be considered within the weakly bending approxima- 
tion for WLCs with a large persistence length L p /L. 

Several generalizations of our approach are possible: 
As alluded to in [16], our model is not in principle re- 
stricted to a pair of cross-linked filaments, but should 
be generalizable to describe the tensile elasticity of a 
stretched, weakly bending WLC bundle, possibly with 
random and/or reversible cross- links. In order to take 
into account non-affine deformation of cross-linked bun- 
dles, we may have to consider also the shearing of cross- 
links. Detailed analysis of bundles exists, due to the 
relevance for actin networks, mostly for reversibly cross- 
linked and extensible, semiflexible polymers [29, 30J. Par- 
ticularly for bundles, the effect of excluded volume inter- 
action, neglected in this work, remains to be explored. 

Apart from activities in this realm, major recent re- 
search efforts are devoted to the impact of structural 
inhomogeneities caused by the (local) breaking of com- 
plimentary base-pair bonds ("unzipping" or denatured 
"bubbles") on the elasticity of double-stranded DNA 
[3T| 132] . A class of semi-microscopic models convenient 
for analyzing the thermal denaturation transition, as well 
as the "bubble" statistics and dynamics |33| . focuses on 
the form of the base-pairing interaction, but does not 
account for the polymers' conformational degrees of free- 
dom, which determine certain DNA properties [33]. In 
the breathing DNA model [35]; two discrete chains (con- 
sisting of interacting "beads" ) with bending and stretch- 
ing rigidity interact via the pairing energy of compli- 
mentary bases, represented by a Morse potential. An- 
other semi-microscopic model, amenable to a transfer 
matrix method, considers a discrete WLC model for the 
chain conformations, coupled to an one-dimensional Ising 
model describing the internal base-pair states |36| . In the 
context of denaturation of DNA, it would be interesting 
to extend our model to reversible cross-linking, in order 
to study the coexistence of "ladders" (cross-linked strand 
sections) and "bubbles" (open sections). A first, obvious 
step towards this direction will be to address with our 
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Appendix A: Structure of the projector sum P 

The matrix representation of the cross-link projector 
sum P, cf. Eqs. (25), has the following structure: 




/ 1 ... 

o '•• 

: 1 







-1 o\ 
• •' : 



-1 







-1 



1 



o .■' 

-1 

V o ... 



'•. 

1 



2Nx2N 



0/ 

(Al) 



Appendix B: Trace of powers of the projector sum 

By decomposing mode indices m 6 {1, 2, . . .} according 
to the block structure into 

rn = 2/iN — r with the definitions (Bl) 

/x := € {1, 2, . . .} (block index), 

r £ {0, . . . , 2N - 1} (index within a block), 

p(r) :— \ \ € {0, 1} (quadrant within a block), 

the entries of P, cf. Eqs. (25 1 and (All), can be encoded in 



product form, the first two factors indicating the location 
of nonzero entries, the last factor the sign, 

Pmimz — ^mi — mi ,2~ZN ^mi+m2,2ZJV (B2) 
= (5 rur2 + 5 rit2N - r2 ){l - 5 ru ZN)(-l) Pl+P2 ■ 
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Then, with the diagonal matrix C, cf. Eq. (18), we write 
the trace of a power k of C~ 1 UU T , cf. Eq. (231, as 



(B3) 



(f)~\r(C-W)* 

oo 

= C ™i (^Wi-m2,2ZiV — <5mi+m 2 ,2ZAr) 1 

mi , . .. ,m£ — 1 

X C ml {3m k -m 1 ,2ZN — ^m fc +mi ,2ZJV ) 
oo 27V-1 

- E E 

fti,...,/i fc =l ri,...,r fc =0 

X C^iV-n^.ra +<Jr 1 ,2iV-r 2 )(l ~ ^.Ziv) (-l) Pl+Pa 
X ^ 2 N-r 2 ( S r2,r 3 + Sr 2> 2N-r 3 )(l-S r2 ,ZN)(-l) p2+P3 
X . . . 

x %Vr k 0Wi + S rk , 2N - ri )(l - 6 rh<ZN )(-iy> + »\ 

Due to the symmetry of the constraints, we are left with 
the summation over one of the rj without and N. If 
we split this sum according to the constraints and to the 
eigenvector structure mentioned above as 



2JV-1 



JV-l 



(B4) 



r=0 



we finally arrive at 



^1 u-K- Vr'V 



(B5) 



oo N-l 



E E ( c 2(/xi-i 



cr 1 



jLll,...,/iH.=l Z = l 
-1 



1 C 2(,u fc -l)iV+i + C 2^ k N-l 



N-l / oo \ ' 

E ( E ( C 2(m-1)JV+J + C 2pJV-i) J 
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Appendix C: Continuous cross-linking 



Here, we present the general result for the limit of 
continuous cross-linking dealt with in Sees. |III B| and 
IIVBI The incremental extension due to this kind of cross- 



linking computed from Eq. ( 34 ) is 

(Ax) 
L 

sing" 
9" 



(CI) 



L 

4L. 



sinh g + 

g+ COth yfjr 



v 



cosh g+ — cos g~ 



as a function of the dimensionless force f r 
the dimensionless parameters 



Eq. (11), and 
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9 



^2^2~gL^/n^J r , 



(C2) 



= j2^/2gLy K + f r 



which apart from f r contain the ratio of total inter-chain 
attraction gl? = 2Nk-QTL 2 /a 2 to bending energy n/L. 
(The function in Eq. (CI) is a continuous, real- valued 
function of f r independently of the sign of the outer 
square root's argument.) 

Comparing the force-extension curves for finite num- 
bers of cross-links to each other and to those for con- 
tinuous cross-linking, we find an approximate collapse 
of all curves with the same total inter-chain attraction 
g a = 2NL 2 /a 2 , as shown in Fig. [6] Except for the 
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FIG. 6. (Color online) Cross-link-induced extension incre- 
ment at total cross-link strength g a = 500 and L p /L = 10. 

case of a single cross-link, the force-extension relation 
is over a large range of forces to high numerical preci- 
sion determined by the product of cross-link number +1 
and strength of a single cross-link only. On the basis of 
Eqs. (34) and (27), this apparent scaling can be traced 
back to the rapid decay of the coefficients c" 1 with mode 



index m. The sum over eigenvectors / in Eq. (27) for 
finite iV is dominated by those with entries at the low- 
est modes, and of the set of modes represented by one 
eigenvector, only the lowest mode gives an appreciable 
contribution. Thereby, the "missing" modes m = ZN 



and, finally, the deviation from the series in Eq. (34) are 
negligible but for very small N. According to the slower 
decay oc m~ 2 of the stretching contribution to the coeffi- 
cients c" 1 , the approximate scaling must break down for 
strong stretching (and small persistence lengths). In this 
regime, however, the additional extension due to cross- 
links is a small quantity anyway. 
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